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Abstract 

The open spin s XXZ model with non-diagonal boundaries is considered. Within 
the algebraic Bethe ansatz framework and in the spirit of earlier works we derive 
suitable reference states. The derivation of the reference state is the crucial point 
in this investigation, and it involves the solution of sets of difference equations. For 
the spin s representation, expressed in terms of difference operators, the pseudo- 
vacuum is identified in terms of g-hypergeometric series. Having specified such 
states we then build the Bethe states and also identify the spectrum of the model 
for generic values of the anisotropy parameter q. 
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1 Introduction 



Much interest has been focused during the past years on integrable models with bound- 
aries [U, |2], further invigorated by advances on problems related to condensed matter 
physics and statistical mechanics, with immediate applications and possible experimental 
realizations (see e.g. [H H]), and also by studies on the understanding of D-branes by 
means of the boundary conformal field theory [51 E] ■ 

The main objective of the present work is the study of the open XXZ model with 
integrable non-diagonal boundaries for the spin s representation. The spin s XXZ model 
[7] with periodic boundary conditions was investigated in [8], whereas studies concerning 
the model with toroidal and generic twisted boundary conditions were presented in [91 [10] • 
Until recently, the formulation of the Bethe ansatz equations for open spin chains with 
non-diagonal boundaries was an open problem. The main difficulty arising when trying 
to obtain the transfer matrix eigenvalues and the corresponding Bethe ansatz equations, 
by standard algebraic or analytical Bethe ansatz techniques, is the lack of an obvious 
reference state 'pseudo-vacuum'. However, the problem was solved rather recently for the 
spin | XXZ chain by two different methods [TT 1 [T2 l [T3]. More precisely, in the approach 
described in |TTl 112] such reference state is not a necessary requirement anymore as long 
as q is a root of unity. On the other hand in the method developed in [13] , which is valid 
for all values of q, suitable local gauge transformations are implemented, along the lines 
described in [TH [To] , rendering the derivation of a suitable reference state possible. 

Note that the spectrum of the spin s XXZ model in the presence of non-diagonal 
boundaries was derived in [IB] using the method of [UJ, however only for q root of unity 
(see also [17] for spin-1 type chains with diagonal boundaries). Moreover, in [18] the 
spectrum, Bethe ansatz equations and Bethe states were derived for the so called cyclic 
representation of U q (sl 2 ), and for the q harmonic oscillator. The generalized spin s XXX 
(q = 1) model with non diagonal boundaries was also investigated within the Bethe ansatz 
framework in [T9], [20] • Here using the algebraic Bethe ansatz framework [21 H3j [151 12] we 
are able to specify the spectrum of the spin s XXZ model with non-diagonal boundaries 
for generic values of q, and more importantly deduce the corresponding eigenstates. It 
should be stressed that the derivation of the corresponding eigenstates is of great signif- 
icance allowing for instance the calculation of exact correlation functions [221 123] . The 
diagonalization process rests on the identification of an appropriate reference state upon 
which all the Bethe states are built. It is also worth remarking that in the case where 
the spin s representation is expressed in terms of difference operators, the reference state 
may be derived by solving sets of difference equations, and it is expressed in terms of 
g-hypergeometric series. 
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2 The open spin s XXZ model 



Before we proceed with the derivation of the spectrum and eigenstates of the spin s XXZ 
chain with non diagonal boundaries we shall first provide a brief review of the model. 
Consider the Lax operator L(\) G C 2 £g> W 9 (s£ 2 ) being a solution of the fundamental 
algebraic relation 

RMi - A 2 ) Li(Ai) L 2 (X 2 ) = L 2 (A 2 ) L x (Ai) ^(Ai - A 2 ), (2.1) 

where R is the spin | XXZ matrix, solution of the Yang-Baxter equation [TH [22] ■ The 
L operator may be written, with the help of the evaluation homomorphism [24] . which 
maps U q (sl 2 ) — > ZY ? (sZ 2 ), as a 2 x 2 matrix 

i(A) = if ^ 9 U-e-V*D ( 4 - r ',B \ 
V ; 2 ^ (g - g-!)C e"VD - e^^U y v y 

A, B, C, D are the generators of W g (s^ 2 ) and they satisfy the well known defining relations: 

AD = DA = I, AC = gCA, AB = g -1 BA, 

r -i A2_f)2 

DC = g _1 CD, DB = gBD C, B = -. (2.3) 

L J q-q- 1 

This algebra admits various realizations, such as the known spin integer or half integer s 
representation. The 2s + 1 dimensional representation may be also realized by difference 
operators acting on the space of polynomials of degree 2s (\2.7\i (see also e.g. |25j). For 
the particular case where the generators admit the spin | representation the L operator 
reduces to the spin | XXZ R matrix. More precisely, A = D _1 i— ► q G \ B i— > a~ , C i— > o~ + , 
with a z,± being the usual 2x2 Pauli matrices. 

(a) Let us recall the spin s representation of W g (s/ 2 ) in more detail. This is a n = 2s + 1 
dimensional representation, and may be expressed in terms of n x n matrices as 

n n— 1 n— 1 

A = <l ak e kk, C = ^2 &k e k k+ i, B = ^C k e k+1 k (2.4) 

k=l k=l k=l 

where we define the matrix elements (eij)ki = Si k Sji and 

n + 1 ~ / q k — q~ k 

a k = — k, C k = J [k)g[n - k] q , [k] q = _ 1 . (2.5) 

Z v q q 



(b) Another realization of the spin s representation, equivalent to the latter one (12.41) . may 
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be given in terms of the Heisenberg algebra — essentially in terms of difference operators. 
Let X and Z be operators satisfying the following commutation relation 

XZ = ? ZX, q = e*. (2.6) 

It is clear that X and Z can be realized in terms of the Heisenberg algebra elements i.e., 

X = ±e p , Z = ±e 4 , and [p, (2.7) 

then the generators A, B, C, D may be expressed as: 

A = q~ s X, D = g-X" 1 , B = ^^(X^ 1 - X), C = ^— (g^X^ 1 - q- 2s X). 

q — q q — q 

(2.8) 

Finally in the 'space' representation |x) e V^ R > of the Heisenberg algebra one defines 

d 

p x) = z/i7rl x ) » <?l x ) = x l x )- ( 2 - 9 ) 
ox 

then it is clear that for any function of z = e x 

X ±x F(z) = F(q ±1 z). (2.10) 

Recall that our objective here is to study the spin s XXZ chain with non diagonal 
boundaries. As is well known to build the model one first needs to construct the open 
transfer matrix, defined as [2] 

t(X) = tr |^ + (A) T(A) K-(X) T(A)| = tr jif + (A) T(A)| where 

T(X) = L Q1 (X) ...L 0N (X), f(X) = L 0N (X)...L 01 (X) (2.11) 

L is given by ( 12.21) and also satisfies L(X) L(—X) = 2 sinh fi(X + is + |) sinh fi(— X + is + |). 
are solutions of the reflection equation [lj: 

R X2 [X X - A 2 ) K X {X X ) R 21 (X 1 + A 2 ) K 2 (X 2 ) = K 2 (X 2 ) i? 12 (A x + A 2 ) K^XJ R 21 (X 1 - X&.12) 

The general solution is a 2 x 2 matrix with entries [261 127] : 

Ku(A) = sinh/i(-A + zO, ^22 (A) = sinh //(A + if) 

ifi 2 (A) = Kg e sinh(2/xA), i^i(A) = Kq- e sinh(2^A). (2.13) 

Here we shall consider K + = K(—X — i; £ + , k + , 9 + ), and K~(X) = K(X; f~ , «T, Using 
the fact that T satisfies the reflection equation one can show the transfer matrix (12.111) 
provides a family of commuting operators [2|: [t(X), t(X')] = 0. 
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3 Bethe states and spectrum 



We can now come to our main aim which is the derivation of the spectrum and Bethe 
states for the generic open transfer matrix (12.1 ip . We shall apply in what follows the 
approach developed in [13] in order to deduce the eigenstates, the spectrum, and the 
Bethe ansatz equations of the spin s XXZ chain with non-diagonal boundaries. Define 
the gauge transformed L— operators 



In $n 



L n (m\X) = M^_ 1 (A)L n (A)M mn (A) = 

S n (m\\) = M-K-X^i-XlM^i-X)^ I £ ), m» = m + X;^3.1) 




k=l 



where the local gauge transformations are defined as (for a more detailed description of 
the method see |13|) 



M n (A) = I x Xn x Vn J , M n (A)= f * n+1 l n - 1 I (3.2) 

x n = — ^ e - i A t (/ 3 +7) e -*M™ ; y n = —i e -w{P-i) e wn ) anc [ ^ are ^ independent C numbers. 
Note that depends on the choice of representation and it will be specified later on. For 
the spin | representation in particular g = 1 [13] (for a relevant discussion see also [18J). 
The transfer matrix (12.111) can be rewritten with the help of the aforementioned gauge 
transformations (13. ip as (see also Appendix and [T3"l ITS]) 



t(A)=tr |K + (A)T(A)j, (3.3) 

where K , T are the gauge transformed matrices (see Appendix). 

To diagonalize the transfer matrix (12.1 ip one first needs a suitable pseudo-vacuum 
state of the general form: 

n<"»> = ^ = M n m) (3-4) 

Wn is the local pseudo-vacuum annihilated by the action of the lower left elements of 
the transformed L matrices (I3.ip . i.e. 

%, % = 0. (3.5) 

From the action of 7„ on the local pseudo-vacuum the following constraint is obtained 

- x mn+1 (e^VA n - e^VDn) + x mn _ 1+1 (e"VD n - e~^q^k n ) 

+e-» X x mn+1 x mn _ 1+1 (q - g- 1 ) C n - e»\q - q' 1 ) B n l = 0, (3.6) 
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a similar constraint arises from the action of 7^ on the local pseudo-vacuum. Notice that 
the transformed non-diagonal elements of (1A.5jl . (1A.7jl acting on the pseudo- vacuum state 
(13 .4p must satisfy: 

K+(m°\X) = K+{m°\\) = , K^(m \X) = (3.7) 

where the integers m° and m are associated to the left and right boundaries respectively. 
Solving the latter constraints (13. 7p one can fix the relations between the left and right 
boundary parameters (see also [TT1 [TBI HE])- Note however that the problem was solved 
in [12] in the more general case with no constraints between left and right boundary 
parameters, however only for q root of unity. 

We may now specify the exact reference state for the spin s XXZ model in both real- 
izations ([22D and (I2TB1 . 

(a) We associate each site of the chain with the spin s n representation and we define 
the local pseudo-vacuum state as 



J 

(m) -^^ (m ' n) // n) (3-8) 



8=1 

fi^ is the j„ = 2s n + 1 column vector with zeroes everywhere apart from the I th position. 

Solving the constraint (I3.6P on the local pseudo-vacuum one can specify the value of 
gn = jn — 1 appearing in the local gauge transformations, and we also acquire recursion 
relations among the coefficients w < f l,n \ which read as: 

w\ m ' n) = ± ^ = q - w\™{ n) (3.9) 

(q-q x ) G_i x mn _ 1+1 

with normalization = 1. The solution of the later recursion formula provides the 

exact form of the coefficients: 
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(m,n) _ ( 1] ^ TT ~ k ^ 



(b) Let us now consider the difference realization of the spin s representation (12.81) . In 
this case the pseudo- vacuum will be expressed as a polynomial of the local spin parameter 
associated to the n th site of the chain 

W M = 0M( z W) = J24 m ' n) (z^) 1 - 1 . (3.11) 
1=1 
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As in the previous case the aim is to determine the coefficients of the polynomial. The 
solution of the constraint ( 13. 6ft this time leads to sets of difference equations for the spin 
variables z^ n \ namely 

(z^x mn+1 q~^ +1 + l) 0("%zW) - (z^x mn _ 1+1 q± + l) ^(g-i^)) = . (3.12) 

It is worth stressing that difference equations of the type (I3.12p occur also in the study 
of the Azbel-Hofstadter problem (see e.g. [281 [22]). Although for our purposes here 
relations (I3.12p are sufficient we derive, for the sake of completeness, the explicit form of 
the polynomial. In any case, the derivation of the pseudo-vacuum is of great significance 
per se leading to the explicit from of the Bethe states as we shall see. In addition we 
wish to stress the connection of these polynomials with the g-hypergeometric series. Let 
us illuminate this point a bit further. Equation ( 13.121) also provides recursion relations 
among w < f 1,n \ whose solution gives the coefficients: 



z-i 



w i 



[in ~ k}c 



M = (^^q-^) 1 - 1 II lA \r jl - ( 3 - 13 ) 

k=i 



It will be instructive at this stage to express the polynomial ( - m - ) in terms of g-hypergeometric 
series (see e.g. [30] and references therein). For that purpose it is convenient to rescale 
the local spin variables as: z < - n - ) = — q % +1 x mn _ 1 +iZ^ n \ 

then it follows form (13.131) that the new coefficients of the polynomial, written in terms 
of T- n \ are given by 



l ~ 2 - r -2j„+2 n 2k 



,-fon) _ TT 1 'I 3 1 



w l 



n • ^ 



q q 

k=0 y * 



Let us also introduce some useful notation. Set q = q 2 , and also define the (/-shifted 
factorials as, 

n— 1 k 

(a;q)n= - aq k ), and (a x , a 2 , ■ ■ ■ , a k ; q) n = JjK; q)n- (3.15) 

fc=0 j=l 

Now we may define the so called g-hypergeometric functions in the following manner 

B +i$ a (ai, a 2 , . . . a s+1 ; b x , . . .b s ; q, z) = > — — z , (3.16) 

^ {b u ...,b s ,q; q) k 

it is thus clear that the polynomial ( 13.131) may be written as a g-hypergeometric series 



M( z ( 



n) ) = £ K—f )k ^ n) ) k =i $ o(f j " +1 iP (n) ). (3.17) 



k=0 
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This series obviously terminates because of the presence of the q ^ n+1 term in the (/-shifted 
factorial appearing in the numerator of (13. 13ft . 

We have specified the exact form of the reference state for both realizations of the 
spin s XXZ model (12. 4p . (12. 8p . To derive the spectrum of the transfer matrix it is also 
necessary to consider the actions of the transformed diagonal elements on the pseudo- 
vacuum. Indeed it is relatively easy to show by simply using the transformations (13.11) . 
and the recursion relations (13.1 3D that the action of the diagonal elements on the pseudo- 
vacua, derived above, takes the general form: 

~g{rn) w {ra) = *m n+ l ~ Vm n -1 ^( A ) ^-1) 
^m„_i+l ?/m n _i-l 

where the values of a n , b n , a' n , b' n for each realization are given below 

(a) a n (X) =g- Jn+1 sinh//(A+j n -), b n (X) = q in ~ x sinh^(A - j n - + i) 
a' n (X) =g j "- 1 sinh / u(A + j n ^), b' n (X) = q~' in+1 sinh ^i(X - ] n % - + i) 

(b) a n (X) = sinh//(A+j n -), b n (X) = sinh //(A - j n - + i) 

<(A) =sinh/i(A + j n ^), b' n (X) = sinh/x(A-j n ^ + 2). (3.19) 

Our aim now is to solve the following eigenvalue problem 

t(X) V = A(A) * where = B m o_ 2 (Ai) . . . B m o_ 2M (A M ) ^ (m) (3.20) 

\I/ is of course the general Bethe ansatz state. 

The transfer matrix eigenvalues may be deduced by virtue of algebraic relations emerg- 
ing from the reflections equation, which omitted here for brevity (see e.g [131 EH])- The 
exchange relations, arising from the reflection equation, are identical to the ones presented 
in [13] . The action of the diagonal elements on the pseudo- vacuum is naturally modified, 
depending explicitly on the choice of the representation. It may be shown that the state 
\1/ is indeed an eigenstate of the transfer matrix if we impose m = m° — 2M, and then it 
follows that mo = Ng + m° — 2M. Let us also define 

f n (A) = a n (X)a' n (X), h n (A) = 6 n (A) b' n (X) 

ts+( \\\ \\\ , sinh^7 + m + 2A + i)sinhz/x ~ + 

K7 (m A) = K7 (m A) H ; : — ^- rKT [m A , 

1V 1 ; sinhi^(l + m + 7)sinh/i(2A + z) 4V 1 ; 

ts+< \\\ sinh^(m + 7)sinh^ ~ + 

KJ (m A) = . - . 7- r — — — -KJ (m A) , (3-21) 
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Kf± are given in the Appendix. Finally the spectrum may be written as 

A (A) = feKiA)irr(m |A)nf n (A) JT ^MA + A,-) sinh//(A - Aj - 
KJ V 1V 1^1^ ' 11 n ^ J-J- sinhu(A + A 7 - + z) sinhu(A-A 

n=l j=l v v 

+ ^K|A)A 4 -(m |A)nhn(A) TT ^hMA + A, + 2j)smh^(A - A^+i) x 

4V i ^ 4 v ui y j_x n V / xi smh /i(A + Aj + smh/i(A-A i ) / v 7 

provided that certain unwanted terms appearing in the eigenvalue expression are vanish- 
ing. This is achieved as long as A,'s satisfy the Bethe ansatz equations, which are written 
below in the familiar form: 

A" 1 + (m°|A i - |) iff (m \\i _ f ) tt f n(Aj - §) ^4 sinh/i(Aj - A, + i) sinh/i(A; + A^ + i) 
Kl(m \Xi - |) AT 4 "(m |Ai - |) £^ h n (A; - |) |Jj sinh/i(A; - A,- - i) sinh/i(Ai + A 3 - - i) 

(3.23) 

if^, f n , h n are defined in (13.211) . and this completes the derivation of the spectrum and 
Bethe states of the spin s XXZ chain with non-diagonal boundaries. 



4 Discussion 

Let us briefly recall the main findings of the present investigation. We have been able 
to derive the spectrum and Bethe ansatz for the open spin s XXZ chain with special 
non-diagonal boundaries and for generic values of the anisotropy parameter. Note that 
the spectrum ( 13.221) . and Bethe ansatz equations ( 13.231) for the spin s XXZ model with 
non-diagonal boundaries are similar with the ones found in [16], although the investigation 
in [16] was restricted for q root of unity. More importantly in the present study we were 
able to specify the Bethe states of the model by means of suitable local gauge transforma- 
tions. The first important step was the derivation of an appropriate reference arising as 
a solution of certain difference equations ( I3.9p . ( I3.12p . Then within the algebraic Bethe 
ansatz framework we built the corresponding Bethe states and derived the spectrum and 
Bethe ansatz equations. 

Acknowledgments: This work was supported by INFN, and the European Network 
'EUCLID'; 'Integrable models and applications: from strings to condensed matter', con- 
tract number HPRN-CT-2002-00325. 
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A Appendix 



The local gauge transformations may be expressed as: 

M n (X)=(x n (X), y„(A)), M„(A)= (x n+1 (X), Fn-i(A)) (A.l) 



and 



X,(A): j 6 ^ ), Y n (X)=( e ^ Vn I. .r„ :,-nr-'>". //„ = //,,r"'" 



(A.2) 



It is also convenient to introduce the matrices 



with 



X n (A) = — ^(-e" A , ar n ), y n (A) = — ^(> A , - y n 
X„(A) = ( - e"\ x n ) , Y n (X) = (> A , - y n ) . (A.4) 

X n - Vn-2 V ' Xn+2 ~ Vn ^ ' 

The transformed K and T matrices are expressed as: 



K + (m\X) 



K+(m\X) K+(m\X) \ Y m (-X)K + (X)X m (X) Y m {-X)K+(X)Y m _ 2 {X) 

K%{m\X) K£(m\X) J \ X m {-X)K + {\)X m+2 {X) X m (-X)K + (X)Y m (X) 

(A.5) 



and 
f(A) 



A m (X) B m {X) \ y m „ 2 (A)T(A)A m (-A) Y m (X)T(X)Y m (—X) 
C m (X) V m (X) J ' \ X m (X)T(X)X m (-X) X m+2 (X)T(X)Y m (-X) 

(A.6) 

Similarly to [13], one defines the transformed K~ matrix as 

K-n\x) = ( ^ (/|A) ) = ( Yi-2(X)K-(X)X l (-X) y(A)^-(A)y(-A) 

^ K 3 -(/|A) KZ(1\\) J \ X l (X)K-(X)X l (-X) X l+2 (X)K-(X)Y l (-X) 

(A.7) 

with I = m + J2n=i9n- 
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